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GLOBAL EXISTENCE AND REGULARITY OF SOLUTIONS FOR 

ACTIVE LIQUID CRYSTALS 


GUI-QIANG GHEN, APALA MAJUMDAR, DEHUA WANG, AND RONGFANG ZHANG 

Abstract. We study the hydrodynamics of active liquid crystals in the Beris-Edwards 
hydrodynamic framework with the Landau-de Gennes Q-tensor order parameter to de¬ 
scribe liquid crystalline ordering. The existence of global weak solutions in two and 
three spatial dimensions is established. In the two-dimensional case, by the Littlewood- 
Paley decomposition, the higher regularity of the weak solutions and the weak-strong 
uniqueness are also obtained. 


1. Introduction 

Liquid crystals are classical examples of mesophases that are intermediate between solids 
and liquids (c/. [H]). They often combine physical properties of both liquids and solids, 
and in general liquid crystals can be divided into thermotropic, lyotropic, and metallotropic 
phases, according to their different optical properties. Nematic liquid crystals are one of 
the most common liquid crystalline phases; nematics are complex liquids with a certain 
degree of long-range orientational order. That is, the constituent molecules are typically 
rod-like or elongated, and these elongated molecules flow about freely as in a conventional 
liquid but, whilst flowing, they tend to align along certain distinguished directions (c/. 

mm)- 

There are several competing mathematical theories for nematic liquid crystals in the lit¬ 
erature, such as the Doi-Onsager theory proposed by Doi m in 1986 and Onsager [H] in 
1949, the Oseen-Frank theory proposed by Oseen |45] in 1933 and Frank |18] in 1958, the 
Ericksen-Leslie theory proposed by Ericksen [15] in 1961 and Leslie |32] in 1968, and the 
Landau-de Gennes theory proposed by Gennes [20| in 1995. The first one is a molecular ki¬ 
netic theory, and the remaining three are continuum macroscopic theories. These theories 
can be derived or related to each other, under some assumptions. For instance, Kuzzu- 
Doi |30| and E-Zhang |14] formally derived the Ericksen-Leslie equation from the Doi- 
Onsager equations by taking small Deborah number limit. Wang-Zhang-Zhang m justi¬ 
fied this formal derivation before the first singular time of the Ericksen-Leslie equations. 
Wang-Zhang-Zhang [62| presented a rigorous derivation of the Ericksen-Leslie equations 
from the Beris-Edwards model in the Landau-de Gennes framework. Ball-Majumdar |3| 
and Ball-Zarnescu [1] studied the differences and the overlap between the Oseen-Frank 
theory and the Landau-de Gennes theory. See [MIIMIEZI for further discussions. 
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Active hydrodynamics describe fluids with active constituent particles that have collec¬ 
tive motion and are constantly maintained out of equilibrium by internal energy sources, 
rather than by the external forces applied to the system. In particular, when the particles 
have elongated shapes, usually the collective motion induces the particles to demonstrate 
orientational ordering at high concentration. Thus, there are natural analogies with ne¬ 
matic liquid crystals. Active hydrodynamics have wide applications and have attracted 
much attention in recent decades. For example, many biophysical systems are classified 
as active nematics, including microtubule bundles [53], cytoskeletal hlaments [29|, actin 
hlaments |6], dense suspensions of microswimmers |63j . bacteria m , catalytic motors [H] , 
and even nonliving analogues such as monolayers of vibrated granular rods |38| • For more 
information and discussions, see |5l[T3l[26l[271ll9l[5T] and the references therein. Active ne¬ 
matic systems are distinguished from their well-studied passive counterparts since the con¬ 
stituent particles are active; that is, it is the energy consumed and dissipated by the active 
particles that drives the system out of equilibrium, rather than the external force applied 
at the boundary of the system, like a shear flow. Consequently, active dynamics are truly 
striking, and many novel effects have been observed in active systems, like the occurrence 
of giant density fluctuations [10lll2l[50], the spontaneous laminar flow [2H|39l[5^, uncon¬ 
ventional rheological properties [I71[25l[55|, low Reynolds number turbulence |271[63]) and 
very different spatial and temporal patterns compared to passive systems [8ll21ll40ll41ll52j 
arising from the interaction of the orientational order and the flow. 

In this paper, we use the Landau-de Gennes Q-tensor description that is one of the most 
comprehensive descriptions, which describes the nematic state by a symmetric traceless 
3x3 matrix, the Q’t^nsor order parameter with five independent degrees of freedom if 
the spatial dimension is three. A nematic phase is said to be (i) isotropic if Q = 0, (ii) 
uniaxial if Q has a pair of degenerate non-zero eigenvalues, and (hi) biaxial if Q has three 
distinct eigenvalues. In particular, a uniaxial phase has a single distinguished direction of 
nematic alignment, and a biaxial phase has a primary and secondary direction of preferred 
alignment. We remark that two-dimensional Q-tensors have been used to successfully 
model severely confined three-dimensional nematic systems that are effectively invariant 
in the third dimension. 

In particular, we consider the following hydrodynamic equations that model spatio- 
temporal pattern formation in incompressible active nematic systems (see [221126] ): 


'dtQ + {u • V)Q + Qn-nQ- X\Q\D = FH, 

< {dt+Ui3di3)Ua + daP - tJ-AUa = di3ral3 + dpaa^, ( 1 . 1 ) 

V • tt = 0, 


where u G d = 2 or 3, is the flow velocity; P is the pressure; Q is the nematic tensor 
order parameter that is a traceless and symmetric dxd matrix; ^ > 0 denotes the viscosity 
coefficient; F“^ > 0 is the rotational viscosity; A G M stands for the nematic alignment 
parameter; D = ^(Vu-I-Vu''') and II = ^(Vu—'S/u^) are the symmetric and antisymmetric 
part of the strain tensor with {Vu)a 0 = djSUa- Hereafter, we use the Einstein summation 
convention, he., the repeated indices are summed over, and a, /3 = 1,2,... ,d. The time 
variable is t > 0, the space variable is x = {xi,... ,Xd), and djs = Moreover, the 
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molecular tensor 

H = KAQ - -(c - c,)Q + b{Q^ - -^^d) - cQtiiQ^) 

describes the relaxational dynamics of the nematic phase and can be obtained from the 
Landau-de Gennes free energy, i.e., Hap = — ^ , where 

^ = / (^(c - c*)tr(Q2) _ ^tr(Q3) + ||tr(Q2)|2 ^ :||vQ|2^d^, 

with K the elastic constant for the one-constant elastic energy density, c the concentration 
of active units and c* the critical concentration for the isotropic-nematic transition, and 
k > 0 and 6 G M are material-dependent constants. We note that the analysis of this 
paper holds for all real 6, but b is usually taken to be positive in the literature. In what 
follows, we set K = k = 1 for simplicity of notation. The stress tensor Gap is 

C^aP = CTap + (^ap, 

with 

(^ap — ^\Q\Hap T QaSHsp HasQsp} ^ aP — ^*C QaPi 

where is the elastic stress tensor due to the nematic elasticity, and is the active 
contribution which describes contractile or extensile stresses exerted by the active particles 
in the direction of the director field (a* > 0 for the contractile case, and cr* < 0 for the 
extensile case). The symmetric additional stress tensor is denoted by 

To/3 ■ Q 0 ^Q)ap dpQ'ySdaQ'yS- 

In the rest of this paper, we consider the case when c > 0 is a constant. We set 

a = i(c-c*), K = a,c\ 

Then system (|l.ll) becomes 


'dtQ + (u • V)Q + Qn-aQ- \\Q\D = TH, 
dtu + (u • V)n + VP - ^lAu = -V • (VQ © VQ) - AV • {\Q\H) 

+V • {QAQ - AQQ) + kV-Q, 


y ■u = 0, 


( 1 . 2 ) 


with 

H = AQ-aQ + b{Q^ - ^^^Id) - cQ tT{Q^), 

the constants c > 0, T > 0, /r > 0, a, 6, A, k G M, and (x, t) G x M+. 

Regarding the related mathematical contributions in the Q-tensor liquid crystal frame¬ 
work, Paicu-Zarnescu [iGlIlTj proved the existence of global weak solutions to the coupled 
incompressible Navier-Stokes and Q-tensor system for d = 2,3, as well as the existence 
of global regular solutions with sufficiently regular initial data for d = 2. Wilkinson [64] 
obtained the existence and regularity of weak solutions on the d-dimensional torus over a 
certain singular potential. In m, Feireisl-Rocca-Schimperna-Zarnescu derived the global- 
in-time weak solutions in the Q-tensor framework, with arbitrary physically relevant initial 
data in case of a singular bulk potential proposed in Ball-Majumdar [3]. Wang-Xu-Yu [60] 
established the existence and long-time dynamics of globally defined weak solutions for 
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the coupled compressible Navier-Stokes and Q-tensor system. See [7] and the references 
therein for more results and discussions. 

We study the active system (II.2p to establish the existence of weak solutions in two and 
three spatial dimensions, along with the existence of regular solutions and the uniqneness 
of weak-strong solutions in the two-dimensional case, motivated by the work of Paicu- 
Zarnescn [461147] for the passive system. Since we are dealing with the active system, 
we need to conquer some new difficulties. Firstly, by using the general energy method, 
we obtain a priori estimates for the system o, based on some crncial cancellations. 
Those cancellations turn out to be very important in the proof of the existence of weak 
solutions in d = 2,, 3, higher regularity and the uniqueness of weak-strong solutions 
in However, due to the appearance of the active term we can only obtain an 

energy inequality, instead of the perfect Lyapunov functional for the smooth solutions 
of the system. Here we mention that the symmetry and traceless properties of the ex¬ 
tensor play a key role in the validity of the cancellations (see also Appendix [^. Also the 
property of the Q-tensor (lA.ip is very important in order to derive the //^-estimate for 
the Q-tensor in Proposition 12.21 since the bulk potential in the Landau-de Gennes energy 
density (the terms independent of VQ) is not always positive. Additionally, in order to 
obtain the weak solutions, we need to add the extra terms —edaQ- VQi 3 .y)\u - VQ\ and 
eV • (Vu|Vttp) to the system to control some non-vanishing terms in the energy estimates 
for the approximate system (13.6p . due to the nonlinearity of the terms: 
and (see §3 for the notations). In ^2]-^ the cancellations for 

these terms work very well. However, in ^ when we seek the regular solutions in the 
cancellation for the terms, X\Q\D and AV • {\Q\H), does not hold perfectly as before. We 
nse the Littlewood-Paley decomposition to reduce these terms to two new terms that can 
be controlled by the previons cancellation idea (see Appendix B for the details). We also 
need to pay more attention to the higher order terms of Q in the elastic stress tensor of 
the system. 

The rest of the paper is organized as follows: In §2, we obtain the dissipation principle 
and a priori estimates. In §3, we establish the existence of weak solutions in d = 2,3. 
In §4, by nsing the Littlewood-Paley decomposition, we restrict ourselves to the two- 
dimensional case and achieve the higher regularity of the corresponding weak solution. In 
§5, we show the nniqueness of the weak and strong solution in with suitable initial data. 
In Appendix A, we provide some important preliminary estimates that we nse extensively 
in this paper. In Appendix B, we provide the detailed estimates for ineqnality (14.61) . 

2. The Dissipation Principle and A Priori Estimates 

In this section, by nsing the energy method, we derive the dissipation principle in 
Proposition 12.11 and the a priori estimates in Proposition 12.21 for system (II.2p . 

For the sake of convenience, we first introduce some notations. We denote H^, with 
k > 1 integer, as the Sobolev space that consists of all functions v in L^(M'^) such that 
is in L^(M'^) for every mnlti-index p = (z^i, • • • , v^), 0 < |p| < k, where D‘' := is 

the distributional derivative. The space is eqnipped with norm || • \\fjk defined by 
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The space H with k > I integer, is defined as the dual spaces of H^, equipped with 


the norm: 




sup \{v,(p)\, 




.<1 


where (•, •) stands for the inner product in L^. For example, if a and b are vector functions, 
then 

(a, 6) = / a{x)-b{x)dx, 

Jm.d. 

and if A and B are matrices, then 

{A,B)= [ A-.Bdx 

with A: B = tr(^i?). We denote by Sq C the space of symmetric traceless Q-tensors 

in d-dimension, that is, 

5^ := {q G = Qpa, tr(Q) = 0, a,/3 = !,■■■ ,d}. 

We define the norm of a matrix by using the Frobenius norm denoted by 

IQI := Vtr(Q2) = sjQapQap- 

With respect to this norm, we can define the Sobolev spaces for the Q-tensors, for example. 


H 


lnn>d cd 


,S^)-.= \q ^ : [ (|Q(x)|2 + |VQ(x)|2)dx < oo 

I JR'i 


We also denote jVQP := dsQa/3dsQai3 and jAQP := 

Let us denote the Landau-de Gennes free energy for the nematic liquid crystals (c/. |20j ) 
by 

F(g) := + f IQP - + l\Q\") dx. (2.1) 

Moreover, by adding the kinetic energy to F((5), we denote the energy of system (ll.2p by 

E{t) := F(Q) + l [ \u\‘^dx. (2.2) 

Proposition 2.1. Let {Q,u) be a smooth solution of system (11.21) such that 

Q G L“(0, T; n T; (2.3) 

and 

u G L°°{0,T]L'^{R‘^))r\L'^{0,T-,H^{R^)) (2.4) 

for d = 2,3. Then, for any given T > 0, we have 

—E{t) + ^ f |Vu|^dx + r [ tic{H‘^)dx < C{k, n) f \Q\‘^dx for any t G {0,T). (2.5) 
dt 2 J^d 

Proof. We take the summation of the first equation in ()1.2I) multiplied hy —El and the 
second equation in ()1.2p multiplied by u, take the trace, and then integrate by parts over 
to hnd 

4 [ (J|VQP + ^|Qp-'^tr(Q3) + ^|Q|4 + l|^|2)^^ + ^||V^x||22+r f tTc{H^)dx 
dt Jwd '2 2 3 4 2 Imd 
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= (u • VQ, AQ) - (u ■ VQ, aQ - b{Q^ - + cQ\Q\^) - (QQ - QQ, AQ) 

a 


+ (ilQ “ clQ — b{Q^ — 


tr(Q2) 


ld) + cQ\Q\^)-X{\Q\D,H) 


- (V ■ (VQ O VQ),u) + X{\Q\H, Vu) + (V ■ (QAQ - AQQ),u) - k(Q, Vu) 

8 

= - k(Q, Vtt) 

i=l 

<^nvuni+c(K,^mil, 

where, in the last inequality, besides Cauchy’s inequality, we have used that X 2 = 0 (since 
V • tt = 0), X 3 + Zg = 0 (by Lemma [A.ip . Zi + Zg = 0, Z 4 = 0, and Z 5 + Z 7 = 0, as shown 
below: 

Zi + Zg = (u • VQ, AQ) - (V • (VQ 0 VQ), u) 

— dadi^QS^'^adx 

= 0 , 

and, by the fact that Q is symmetric and 0 is skew symmetric, 

Z 4 = (OQ - Qn, aQ - b{Q^ - + cQIQp) 


d 

= -{QQ + Qn, aQ - b{Q^ - ^-^^Id) + cQ\Q\^) 

+ 2(QQ, aQ - 6(Q2 - + cQIQp) 

= 0 , 


and 


Z 5 + Z 7 = A(|Q|F, Vn) - X{\Q\D, H) = X{\Q\H, Vu) - X{\Q\H, D) 
= A(|Q|F, Vu-D) = X{\Q\H, n) = 0. 


□ 


Remark 2.1. For the passive system considered in [36], a perfect Lyapunov functional is 
available. However, for the active system as analyzed here, only an energy inequality (12.Sp 
is obtained above, which is not a Lyapunov functional in general. 

Based on Proposition 12.II and Gronwall’s inequality fLemma lA.dp . we have the following 
a priori estimates. 

Proposition 2.2. Let (Q, u) be a smooth solution of system (11.21) in R'’*, d = 2,3, with 
smooth initial data (Q(x),u(x)). If {Q,u) G x L?, then, for any t > 0, 

||Q(t,-)ll//i <Cie^^*(||Q||^i + ||h||i,), 


( 2 . 6 ) 
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and 

2 Olliz + ^ J •)lli2iis < (73(11(511^1 + ||«|li2)e‘^2i _|_ (7^^ ^2.7) 

where constants Ci, 1 < i < 4, depend on (a, b, c, k, /i, A, T,Q,u). 

Proof. Prom the energy estimate in Proposition 12.11 and the symmetric property of Q, we 
have 

j^E{t) + f^WVuWl, + r\\AQ\\l, + a^rWQWl, + c^r\\Q\\% 

+ f \ d)'^)dx 

jRd d 

< C(k, p) WQWl^ + 2ar(AQ, Q) - 2acr||Q||^4 - 26r(AQ, Q^) ( 2 . 8 ) 

+ 26cr(Q tr(Q2), q 2) + 2a6r(Q, Q^) + 2cr(AQ, Q tr(Q2)) 

4 

= C{k,p)\\Q\\1, - 2ar||VQ||i. - 2acr||Q||i4 + Y.Z,. 

i=\ 

We now derive the estimates for Xj, 1 < i < 4. First, we have 

Zi = -26r(AQ,Q2) < ^\\AQ\\l,+C{h\T)\\QtL,. (2.9) 

From (jA.ip . we have the following estimates for X 2 and X 3 by choosing an appropriate 
e > 0: 

X 2 = 26cF(Qtr(g2)^Q2) ^ ^bcT [ tT{Q^)\Q\‘^dx 

<2|6|cF [ (£|g|4 + i|Q| 2 )|Q| 2 rf^ (2.10) 

jRd 4 e 

= A||(}||6. + C(6lr)||Q||l., 

and 

l 3 = 2abT{Q,Q^) = 2abr [ tr(Q3)dx < (7(a, 6 ,F)(HgH^^ + ||Q||| 4 ). (2.11) 

jR<i 

Moreover, we observe that 


I4 = 2cF(Ag, gtr(g^)) = 2 cF f d^^QajiQaii^^{Q^)dx 

Jr<^ 


= -2cF / d...^Qapd..^Qap^r{Q^)dx - 2cF / d...^QaiiQalid..f'^^{Q^)dx 

JRd jRd 

=-2cF / |Vg|2|gpdx-cF / |Vtr(g2)pdx < 0. 

JR'* JRd 

Combining (j2.8p with (j2.9l) — (j2.12p . we have 


jE{t) + |||Vrx||i. + ^||AQ||2, + ^IIQIlie + b^T - ^-^l,)^)dx 

< C{a^b^c,K,p,T){\\VQ\\l, + Iigiii^ + ||g||i4). 


( 2 . 12 ) 


(2.13) 
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Here we should clarify that, since the potential terms, he., the Q-terms without deriva¬ 
tives in E{t) do not always sum to a positive quantity, (|2.13l) does not always yield the 
desired estimates for Q directly. However, we can deal with this issue as follows. 


Case I: a > 0 and sufficiently large. By the property of Q in (lA.ip . we have 



3 V 4 


+ Tl^r > - o 7ltr(Q')r + -tr(Q^) +^\ q \ 




IQP + 


12 ^ 




(2.14) 


When e > 0 is sufficiently small and a > 0 is sufficiently large, both f — ^ and | — 
can be positive. As a result, we can obtain the //^-estimates of Q directly from (12.131) by 
Gronwall’s inequality fLemma IA.4p . 

Case II: For all other a. In this case, the sum of the Q-terms without derivatives in 
E{t) may be negative. Thus, we have to deal with the L^-estimates of Q separately to 
obtain the Lf^-estimates for Q. In fact, we multiply the first equation in (jl.2|) by Q, take 
the trace, and integrate over by parts to obtain 

= - mQWh - - cr||Q|li4 + br [ tiiQ^)dx + X{\Q\D, Q) 

I at jjjd 

< - rllVQIli^ - aVWQWl, - cTWQWl, + C{a^ 6 ^, r)(||Q||2, + \\Q\\l,) 

+ e|A|||Vu||i.+C'(e)|A|||g||i4 

< - rllVQIli^ + elAlllVulli^ + C{\\Q\\l, + \\Q\\l,), 


where e > 0 will be decided later, and C = C{a, b, c, A, T, e). Motivated by Case I, we notice 
that, for any Q G Sq, there exists a positive, sufficiently large constant M = M{a,b,c) 
such that 

0 < ylet + < (V + f)IOt - hr(e=) + jlOI*. (2.16) 

Multiplying (|2.15p by 2M, adding it to (|2.13l) . and letting e = we have 

— [Eft) + M||g||^2) + ^||Vu||^2 + —||AQ||^2 + -^IIQIlie ^2 

<C[\\Q\\l2 + \\VQ\\l, + \\Q\\l,), 
where C = C[a,b,c, k, fj,, X,T, M). 

Then the desired estimates for {Q,u) follow from Gronwall’s inequality 

f Lemma IA.4p . □ 


3. Weak Solutions 

In this section, we prove the existence of weak solutions for system (II.2p with suitable 
initial data for d = 2,3. 
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Definition 3.1. {Q,u) is called a weak solution of system (11.21) with the initial data: 

Q{0,x) = Q{x) e u{0,x) =u{x) € L‘^{R^), V • h(x) = 0 in p'(]R'^), (3.1) 

if {Q, u) satisfies the following: 

(i) QGL^^(R+,H^)nLl^{R+;H^) and G L-(M+; L^) n Lf^,(]R+; i/i); 

(ii) For every compactly supported G C'°°([0, oo) x S'q) and ip G C'°°([0,oo) x 

M'^) with V ■ ip = 0, 

poo p 

/ j {Q \ dtip + FAQ + Q {u ■ VaiVj) — (QH — — A|(5|D) : (fjdxdt 

Jo JR'* 

[ [ (aQ — biQ^ -+cQtr(Q^)) -.(pdxdt— [ Q{x) : ip{0,x) dx, 

Jo d ' / Jmd 

(3.2) 


= F 


and 

roo 


poo p p 

/ / i — u-dtip — u-{u-SIx'ip)+iiVu\V'ip^)dxdt— / h(x) ■-0(0, x) dx 

Jo JmA JRd 

= J J (VQO VQ + A|Q|id-(QAQ-AQQ)- kQ) : VV’dxdt. 


(3.3) 


Theorem 3.1. There exists a weak solution {Q,u) of system (11.21) subject to the initial 
conditions ra, for d = 2,3, satisfying 

Q G Li^,(M+; H^) n LL(M+; H^), u G L^,iR+, L^) n LL(K+; (3-4) 


Before proving Theorem 13.11 we introduce some useful notations: 

(i) is the convolution operator with kernel where y G is a radial 

positive function such that 



= 1 . 


(ii) The mollifying operator n = 1, 2, • • • , is defined by 


R{Jnfm :=l[2-,2"](lel)-^(/)(0, 
where T is the Fourier transform. 

(iii) V is the Leray projector onto divergence-free vector fields, i.e., 

V : L‘^{n) ^ Id = |w G : V • w = o| , 

which can be explicitly described in the Fourier domain by the tensor as 

Ww))(0 = (id - ^).^(w)(0. 

Next we prove Theorem 13.11 in the following three subsections. In ^3.11 we construct 
regularized approximate solutions {Q'f'\u'f) to the approximation system (|3.5p . In 113.21 
similarly to Proposition 12.21 we obtain some a priori bounds in ()3.7p . which allow us to 
obtain the convergence result (I3.10p for by the Aubin-Lions compactness lemma. 

Moreover, we can pass to the limit as n goes to infinity to achieve the weak solution (Qs, u^) 
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of the modified system (I3.11|) . In ^13.31 by studying the e —0 limit of the modified system 
(13.111) . we obtain the weak solution to system (11.21) . However, we cannot use the previous 
a priori bounds in (13.7h . since those bounds are not uniform with respect to e. Instead, we 
need to repeat a similar procedure in order to obtain the uniform bounds in (I3.14l) - (l3.15p . 


3.1. Regularized Approximation System. Let us consider the following approxima¬ 
tion system for the active hydrodynamic system (m, followed by the classical Friedrichs’ 
scheme, for any fixed e > 0 and n > 0 (from now on, solution is denoted by 

, m”) for simplicity of notation when no confusion arises): 

+ JniiVJnReU^ ' V) - Jn{VJn{Re^^)JnQ^"^) 

+ JniJnQ^^^VJn{Ren^)) - A J„(| |iP = FJ^R^, 

dtU^ + VJniiVJnReU^ ' V)VJnU^) - fiAVJnU^ 

^ =-eVJnRe{daJnQf^{ReJnU^-'^JnQf^)\ReJnU^-VJnQ^^^\) ( 35 ) 

+eVS/ ■ JnRe{S/JnReU^\'^JnReU^\^) 

-vv ■ o - xvv • 

+VV ■ JnReiJnQ^^^^JnQ^"^ - JnQ^"^) + kVV • JnReQ^"\ 

SQ^^\u’-)\t=0 = {JnReQ.JnReU), 


where 

rW = - ag(") + b{{JnQ^^'>f - i^) _ cj„gWtr((j„gW)2), 

This approximate system can be regarded as a system of ordinary differential equations 
in L^. By checking the conditions of the Cauchy-Lipschitz theorem m, we know that it 
admits a unique maximal solution (g(”),ri”) e C^([ 0 ,T^);x on 

some time interval [0, T„). By simple calculation, {VJn)‘^ ='PJn and = Jn, so that pair 
{JnQ^^\VJnu'^) is also a solution of (j3.5p . By uniqueness, {JnQ^"'\VJnU^) = {Q^'^\vP'). 
Therefore, (g^""^ , li”) also satisfies the following system: 


+ J„(R,«-VgW) - Jn{ReXl^Q^^^ - g(")R.O") - AJ„(|g(")|R,L»-) = 

dtU^ + VJniReU^'^U^) - 

= -eVJnRe{daQ^p^{ReU^ ' VQ^;^)\ReU^ ■ VgW|) 

< +ePV-JnRe(VRe«”|VReU”P) 

-vv • JnRe{yQ^^'> o vgW) _ xvv ■ j„R,(|gW|j„RW) 

+VV ■ j„R,(gWAgW - AgWgW) + kVV ■ {JnRsQ^^^), 
(gW,«-)|i=o = {JnReQ,JnReU), 

where G C^([ 0 , T^); and 




(3.6) 


= Ag^*^) - ag^*") + 6 ((g ("))2 - 


tr((g("’))^ 


^irf)-cg(")tr((gW) 2 ). 


Remark 3.1. It is easy to see that, if g^”) is a solution to system (13.51) . so is (g("-))'''. Hence, 

g(n) ^ ^g(n))T 

in [0,T„] X by the uniqueness of the solution. Moreover, from 
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now on, we will work with the solution of the system (I3.6p with this symmetry property, 
instead of the solution to the system (|3.5I) . 

3.2. Compactness and Convergence as n —> oo for System (|3.6|) . First, we need to 
derive some a priori estimates for the system (13.6p in the following proposition. 

Proposition 3.1. The solution of the system (13.61) satisfies the following esti¬ 

mates, which are independent ofn, for any T < oo: 

sup \\ReU"' ■ VQ^"'^||i3(o^T;L3) + SUp \\^ReU"'\\L*{p,T-,Lfi < C*, 

n n 

sup 11^2(0,T;Ll2)nLoo(0,T;LIlnL4) + SUp || |^) 11^2 (o,T;L 2 ) < C, (3.7) 

n n 

sup ||^^”||l,°°(0,T;L2)nL2(0,T;Hl) < C, 

n 

provided the initial data {Q,u) € x L^. Moreover, if Q ^ Sq, then G Sq. 

Proof. Similarly to the proof of Proposition 12.11 we sum up the hrst equation in (13.61) 
multiplied by — — c(5("'Hr(((5*'"'^)^)) and the second equation 

multiplied by n"', take the trace, and integrate by parts over to obtain 

Jt L f “ ^tr((QW)3) ^ £|g(n)|4 ^ 1 |^n|2^ 

+ fi\\Vu^l2 +r\\AQ^^^\\l2 + a^T\\Q^^^\\l2 + 2acr\\Q^^^\\l4 
+ 6"r||Jn(Q(’^^)'|li2 + c2r||j„(QW|Q(-)|2)||2^ + 2ar||vg(’^)|li2 
+ £\\ReU^ ■ + eWVRsum^ 

- + (3.8) 

+ {TniRe^^Q^^^ - -6j„(gW)2 + cj„(gW|g(")|2)) 

- A(|g(-)|ii,ii”, Jn(AgW - aQ(”) + 6(g("))2 - cQ(")|g(’^)|2)) 

+ 2r(-6Ag('^) + a6g(") + 5cJ„(g(")|g(’^)|2), j„(gW)2) 

- (V- (vg(")0vg(”)),i?,u”) + A(|g(’")|Jn^^”\vi2eu”) 

+ (V • (g(")Ag(”) - Ag(”)g(")),i?,n") - K(g("), V7?,u”) 

+ ^(tr((g("))2)irf, j„(Ag(-) -ag(") + 6(g(-))2 -cg(’^)|gW|2)) 

12 

= I1^- 

i=l 

By using cancellations analogous to §2, we have 

2^1 + 218 = 0, X3 + ZiQ = 0, 06+219 = 0, 04 < 0. 

The remaining terms can be estimated as follows: 

02 = (j„(i?,n’^vg(")),6j„(gW)2 - cj„(g(")|g(")|2)) 
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< (^ + + ^11 )iii2 + c(62)||qW||4, 

< ^-\\R,u^ • Vg(")||i3 + C{b\e,T)\\ReU^ • Vg(")||^i 

r^V 

+ ^liA(ei">ig'”'i")iii. + c(62)il(3<”i|li. 

< . vo('‘>iii, + c{b\E,T) (ll•<”lli. + iivo>’‘>iii,) 

+ —liA(gi">iQ<”'i")iii. + c(6")ilQ<”'ili4. 

A = - g‘"*fl.f!"), -6j„(g'">)2 + cj„(g'">|g<“>p)) 

< (f + + ^liA(g<’‘>ig'"f )ii|, + c(62)iig(")|ii. 

<|tivftu”||‘. + c(62,r,g||g<“>||i4 + ^|^||j„(g''‘>|g<”f)lli., 

Xr = 2r( - 6AQ(^) + + 6cJ„(Q(’^)|Q(^)p), J„(Q("))2) 

< ^||Ag<")|||, + A||j„(q(«)|q(»)| 2)||2, +C(a^62_r)(||g(")|ii4 + ||g‘”>|il,). 

Ill = -K(g'’‘>,vi;4u”) < ^l|VR,u”||J. + c(K^^)||g(”)||J, 

<fl|v«”lii. + c(A,M)lig‘”'lli4, 

A 2 = ^(tr((g<"i)2)ij, j„(Ag("> - aQ<"> + 6(g(">)2 - a3(")|g<"i|2)) 
a 

< ^l|Agl”>|li. + A||j„{0(«)|q(«)| 2)||2, +C||g(”)||2, + C||Q<”l||l4. 

Substituting all the above estimates into (13.8|) . we have 
d 3/; r r^r 

-E”(t) + ^i|VK”iii4 + -iiAg(»)|ii4 + —I|j„(gi”>ig'”'i")iii4 

+ . Vg<">||i, + |l|VR4«"||t4 (3.9) 

< c (||g<”i|ii4 + iiQ<”>iit4 + iiu’iii. + livggiii",), 

where C depends on a, b, c, k, F, /i, and s, and 

E^{t) = Q|vg(-)|2 + ||QW|2 - ^tr((QW)3) + ||g(n)|4 ^ 1 |^n|2^ 

Again, by the same reasoning as in Proposition 12.21 E^{t) may be negative. We also 
need to estimate the L^-norm of the Q-tensor separately in order to obtain the desired 
if ^-estimates for Q. Multiplying the first equation in system (j3.6p by Q^"‘\ taking the 
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trace, integrating over by parts, multiplying the result by 2M {M > 0 is sufficiently 
large), and adding it to (13.9p . we have 

I + MlIQ^Ili.) + |||Vu"||i. + ^IIAQ^Ili, 

+ + IWReU^ ■ VQW||33 + £||Vi2,unii4 

< C + ||g(")||i 4 + \K\\l, + ||vg(-)||i.), 

where C = C{a, b, c, n, A, F, e, fj,, M) is a constant, independent of n. 

From the above estimate, along with Gronwall’s inequality iLemma IA.4p . we can con¬ 
clude the a priori bounds in (|3.7p . which are independent of n, for any T < oo. 

In order to prove that g^”^ G Sq, besides the symmetry property mentioned in Remark 
EH it remains to show tr(g("')) = 0. We take the trace on both sides of the first equation 
in system (|3.6p and use that g^”^ = (g^”'^)"'", (0"')''' = —0"', and tr(iA”) = div(u”) = 0 to 
obtain the following initial value problem: 

attr(g(")) + Jn{ReU^ • Vtr(g(”))) = FJ„(Atr(g(-)) - atr(gW) - ctr(g(’^))tr((g("))2)), 
tr(g(”))|t=o = JnRetr{Q) = 0. 

Multiplying the above equation by trpg^”)), integrating by parts over and using 
= g^”^ and the uniform bounds of g^”^ in (j3.7p . we have 

:^l|tr(g("))||i 2 + F||Vtr(gW)||i. = -aF||tr(gW)||i3 - cF / |tr(gW)|2|gW 
at j^d 

< -aF||tr(g("))||i, +C||g(")||ie||tr(gW)||^a||tr(gW)||i. 

< -aF||tr(gW)||i. +C||g(")||^i||Vtr(g("))|||,||tr(gW)||^ 
<^||Vtr(gW)||2,+C7||tr(g(-))||i,, 


thus. 


^l|tr(Q'"l)|||, <C||tr(Q(">)l|i., 

where we have used the Sobolev imbedding, the Gagliardo-Nirenberg interpolation in¬ 
equality in Lemma lA.31 and the Cauchy inequality. Hence, we conclude that trpg*^”)) = 0 
by the initial condition. □ 


Then we can conclude from the uniform estimates in (j3.7p that T„ = oo. In addition, by 
using system (j3.6p and the above estimates, we can compute the bounds for dt{Q^'^\u'^) in 
some L^(0,T; H~^) for large enough N. Then, by the classical Aubin-Lions compactness 
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lemma fLemma lA.2l) . we conclude that, subject to a subsequence, 

^ Q in L\0,T-Hf-J) ioi any 5 e {0,2 + N), 

Q{t) in for anyt G M+, 

gW ^QmLPiO,T-,H^), ^ Q in LP{0,T-HI-J) for any 5 G (0,1 + iV), pG [2,oo], 
inL^{0,T-,H^), ^ u in l 2(0, T;for any <5 G (0,1 + iV), 

^ u{t) in L? for any t G M+. 

(3.10) 

As a result, we can pass to the limit as n goes to infinity to obtain a weak solution 
{Qs,Ue) of the following modified system (for simplicity, we denote {Qe,Us) by {Q,u) 
when no confusion arises): 

'dtQ + Reu • VQ + QR^n - Re^Q - \\Q\ReD = TH, 
dtu + V{R£U ■ Vu) — fiAu 

< = —e'PRe{daQi3'y{ReU ■ VQj3^)\RsU ■ VQI) + eW ■ Rs{VR^ulVRsu\'^) 

-VV ■ Re{VQ O VQ) - AV • VRe{\Q\H) + VV ■ R^iQAQ - AQQ) + kVV ■ ReQ, 
{Q,u)\t=o = {R£Q,Rs'ii), 

(3.11) 

such that 

Q, G H^) n lL(^+; u£ g l£',(m+; l^) n lL(^+; (3-12) 

Remark 3.2. This modified system is obtained by mollifying the coefficients of the Q-tensor 
equation and the forcing terms of the velocity equation, and by adding the extra terms 
given by —£daQp^{u ■ VQp-^)\u ■ VQ| and eV • (VujVup) to the velocity equation. These 
two terms are needed to estimate some bad terms which do not disappear in the energy 
estimates. Moreover, we would like to mention that the above procedure for obtaining 
the solution to system (13.lip follows from the classical Friedrichs’ scheme. We also point 
out that the solutions to (|3.11l) are smooth, because we can bootstrap the regularity 
improvement provided by the linear heat equation to obtain the smooth regularity of Q, 
and bootstrap the regularity improvement provided by the linear advection equation to 
obtain the smooth regularity of u. 

3.3. Compactness and Convergence as e —> 0 for System (13. lip . In this subsection, 
we show that, by passing to the e —>■ 0 limit in system (13.111) . we can obtain a weak solution 
of the system (|1.2p . In order to do so, we need to achieve some uniform bounds for solution 
(Qe,tt£;). Although, as a limit of (Qe^^w”), (Q£,Ue) still satisfies the a priori bounds in 
(1221), we cannot apply these bounds in this step because they are not uniform with respect 
to £. Therefore, we need to hnd new a priori bounds for the system (13.lib . For simplicity 
of notation, we denote {Q^jUs) by {Q,u) when no confusion arises below. 

Applying the same procedure as in lj2l that is, multiplying the first equation in (13.lib by 
—H + 2MQ (with M a sufficiently large positive constant), taking the trace, integrating 
by parts over and adding this to the second equation in ()3.1ip multiplied by u and 
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integrated by parts over with the analogous cancellations as in ^ we have 

Tt X. (2 + 2 

+ ^II^^IIl 2 + -^II^QIlia +-^IIQIlie + ell-ReU • VQII^s + e||Vi?£tt|||4 ^ ^ 

< C{\\Q\\l, + \\VQ\\l, + \\Q\\U, 

where C = C{a,b,c, k, fi, X,T, M) is independent of e. Then, by Gronwall’s inequality 
fLemma lA.4h . we have the following a priori bounds, independent of e, such that, for any 
T < 00 , 

sup ||Q£:||L2(o,T;Ll2)nL°o(0,T;LIlnL4)nL6(0,r;L6) + SUp \\Ue\\L°^(^o,T-,L'^)nL^{0,T-,H^) < C. (3.14) 

e e 

Moreover, for any e > 0, we have 

elliJeW • V(5|||^3 (o,T;L 3) + ^11 Vi2eU||^4(-Q <C (3.15) 

with constant C independent of s. In addition, since {Q^,Ue) satisfies system (|3.11l) . along 
with (I3.14p - (|3.15l) . we can obtain bounds for dt{Qe,Ue) in L^{0,T; In order to do 

this, we need to estimate the H~^-norm. of each of the other terms in system (j3.Iip . aside 
from {dtQe, dfUs). Since the estimates are similar to each other, we just show some tricky 
ones in the following: 

\\X\Q\ReD\\H -2 = sup {X\Q\R,D,ip)<C\\X7u\\L2\\Q\\L2ML^ 

< C\\Vuh2\\Qh2\MH2 < C\\X7uh2\\Qh2, 

WeW ■ R£{'VR£u\VR^u\‘^)\\jj -2 = sup {eW ■ R^{'VR£u\VR^u\‘^),'4!) 

IIV'll;f2<l,div(V’)=0 

0 

= sup {£X/RsU\VRsU\‘^,RsV'4!) < C£\\VRsU\\]^2\\VRsU\\\4 

11^11^2 <l,div(V')=0 
"o 

< C'e||Vu||i2 ||Vi2ert|||4, 

\\XV-VReilQlH^H-^ = sup {XV ■VRe{\Q\H),^l;) 

llbll^ 2 <l,div(i/i )=0 

"o 

sup {X\Q\{AQ-aQ + b{Q^-^-^^ld)-cQtv{Q^)),ReV^I;) 

llbll;f2<l,div(V')=0 “ 

"o 

< C\\Q\\l^{\\AQ\\i2 + \\Q\\l2 + \\Q\\‘j^4 + llQIlie), 

and 

\\kVV ■ ReQWn-^ = sup {kTV ■ R,Q,^Ij) < C\\Q\\l 2 . 

IIV’llif2<l,div(V')=0 
0 

By the Sobolev interpolation inequality fLemma IA.3I) . one has 

||VQ||i3<C||Zl2Q|||^||Q|||^. 


(3.16) 
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From the above inequality, we have 

\\eVRe{dcQp^{ReU ■ VQpr^)\ReU ■ VQ\)\\h-2 

= sup {£T’Re{daQj3'y{ReU-VQi3^)\ReU-VQ\),'lp) 

ll'!/'llH^2<l,div(i/i)=0 

"o 

= sup [edaQji^iReU ■ VQp^)\RsU ■ \/Q\,Rsip) 

IIV'lljj-2<l,div(i/i)=0 

"o 

<Ce\\VQ\\LA\ReU-VQ\\%\\nL^ 

< Ce\\D^Q\\lA\Q\\lA\ReU • VQIliallV’llL- 

< Ce{\\ReU ■ VQIlis + \\D‘^Q\\l 2 + \\Q\\le). 

From all the above estimates, together with the uniform bounds (l3.14jl - (j.S.15jl . we can 
conclude that {dtUe,dtQe) G L^{0,T; H~^). Then, by the classical Aubin-Lions compact¬ 
ness lemma (Lemma IA.2p . we know that there exist Q G H^) n R'^) 

u G L^) n R^) such that, subject to a subsequence, we have 

inL^(0,T;R^), Q, ^ Q in L\0,T; hI-J) for any <5 G (0,4), 

Qe{t) Q{t) in for any t G M+, 

Qe^Q mLP{0,T-,HA, Qe^Q in LP{{),T; H]-^) for any 5 G (0,3),p G [2, oo). 
Usu in L^(0, T; and Us ^ u in L‘^{0,T; for any <5 G (0, 3), 

Us{t) u{t) in for any t G M+. 

(3.17) 

With the above result, we can pass to the limit in the weak solution {Qs,Us) of system 
(j3.1ip . as e —>■ 0, to obtain a weak solution (Q, u) of system (jl.2l) satisfying (j3.2l) - (j3.3l) . In 
the following, we focus on some terms that are not so easy to deal with. 

First, we observe that 

dfsi.Qe^^sRsidjj'lpa^ dfsQ-ysdfi'Aci 

= dii{Qs)^5{Re{dpAa) - dpAa) + {dp{Qs)^S - di3Q^s)d0Aa (3.18) 

converges to zero strongly in L^(0, T;L^). This is owing to the fact that converges 
to zero strongly in L^(0, r;L^), since Rsidjsipa) — djsipa converges strongly to zero in any 
LP{0,T; L’^) (ip is compactly supported and smooth) and Qs is uniformly bounded in 
L°°(0, T; iF^), and the fact that X| converges to zero strongly in L^(0,T;L^), since Qs 
converges to Q strongly in L‘^{0,T; and A is compactly supported and smooth. 

Combining the above facts with the weak convergence of Qs in L^{0,T; H^) in (13.171) yields 

roo r POO p 

I I ^e{9a{Q£)'ysdi3(^Q^')'yS^di3'ijjadxdt — / / do^(^Q^')^sdis{Q£)'ysdi£di3'ijj(ydxdt 

Jo Jr^ Jo Jr^ 
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Moreover, from the strong convergence of in for p G [l,oo) in ()3.17p 

and the uniform bound of Qe in L®(0,T;L®) in (I3.14p . we have 

\Qe?Qe ^ \Q?Q, \Qe\Ql ^ in L\Q,T-,Ll,). 

As a result, we have the following convergence: 

poo p poo p 

X / Re{\Qe\H^) :S/'iljdxdt = X / \Q,\H^ : ReVipdxdt 

Jo JR<i Jo Jr'^ 

= x y‘'"y‘jQ,|(AQ,-aQ, + 6(Q2 + M^£lil^) _cQ,tr((Q,)2)) :R,V^dxdt 

(3.19) 

poo p 

^X / \Q\H-.Vijjdxdt. 

Jo Jr‘^ 

Finally, from the uniform bounds of e||i?£n • in (|3.15l) . along with (I3.16p . we have 

poo p 

e / Re{\ReU-X/Q\{RsU-X/Qi3^)daQi3'y)'J’adxdt 

Jo 4Rd 

: / / {^\RsU ■ X/Q\{RsU ■ VQpj)daQpry)Ri;'ijjadxdt 

Jo JRd- 

:f \\ReU-VQ\\l3\\VQ\\L3\\ReJj\\L--dt (- 320 ) 

0 

T 

Ce [ \\R,u-X/Q\\l3\\D^Q\\l,\\Q\\l,dt 
Jo 


= e 


< e 


< 


< Ce\\ReU ■ VQ |||3(o,T;L3) IIQII L 6 ( 0 ,T;L 6 ) l^{o,t-,h^) 

0 . 


Similarly to the above estimate, by using the uniform bound of s\\X/Reu \\^4 in (13.151) . we 
also conclude that e f (Vi?£M|Vi?£Up) : V'lpdxdt ^ 0 as e —>■ 0. 

Then the proof of Theorem 13.II is completed. 

We remark that, in the proof above, because of the active terms in our system (ll.2jl . 
in order to obtain the convergence of the approximate solutions, we need to establish a 
higher integrability of Q in time and use the Sobolev interpolation inequalities to achieve 
the uniform F7“^-estimates for the extra terms added to the approximate system, which 
is different from the passive case treated in [46]. 


4. Higher Regularity In Two Dimensions 

In this section, we prove that, in the two-dimensional case, system dOI) has solutions 
with higher regularity, provided with sufficiently regular initial data. The result is stated 
in Theorem liTl We mention that we use the Littlewood-Paley decomposition to help us 
improve the higher regularity of the solution. Of course, we can also obtain the higher 
regularity by differentiating the equations k > 1 times in system (|1.2p . However, this 
requires the initial data (Q,u) to be at least in x rather than in x for 

s > 0 in the Littlewood-Paley method. 
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Remark 4.1. For any traceless, symmetric, 2x2 matrix Q, 



which means that H reduces to a simpler form 


H = AQ — aQ — cQtic{Q'^) for x G M^. 


Theorem 4.1. Let s > 0 and {Q,u) € x There exists a global solution 

{Q{t,x),u{t,x)) to system (|1.2I) with the initial conditions: 


Q{0,x) = Q{x), u{0,x)=u{x) 


(4.1) 


such that 


Q G nL-(M+;17^+1 (r2)) 


and 



(4.2) 


where constant C depends on Q,u,a,b,c, and F. 

To prove this theorem, we restrict ourselves to system (11.21) in two spatial dimensions 
and give the a priori estimates for the smooth solutions of this system. Of course, the 
same estimates independent of e can also be obtained, if we use the modified system 
(jd.llj) . whose solutions are smooth as we mentioned in Remark [321 Moreover, we use the 
Littlewood-Paley decomposition to obtain the a priori estimates of the solution. Firstly, 
we apply (see Appendix A for the notations), with g G N, to the equations in system 
o to obtain the estimates of high frequencies of the solution. Secondly, by applying Sq 
to this system, we obtain the estimates of low frequencies. Finally, we achieve the high 
regularity of the solution by combining the high and low frequencies together and by using 
Gronwall’s inequalities fLemma IA.4I1 . 

Proof. We begin with the estimates of high frequencies. Applying Aq to the first equation 
in (m, using Bony’s paraproduct decomposition, multiplying the equation by —AAgQ, 
taking the trace, and integrating by parts over we have 







(4.3) 


q'>q-5 
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+ ^ ^ j^2^'iP^q'Qa-y)j ^^qQa/s) 

q'>q-5 

- X{Aq{\Q\D), AAqQ) + TaiAqQ, AAqQ) + Tc{Aq{Qix{Q^)), AAgQ) 

10 

= E^- 

1=1 

Similarly, we can apply Ag to the second equation in (11.21) . use Bony’s paraproduct de¬ 
composition again, multiply the equation by AgU, and integrate by parts over to obtain 

+ /r||VA,tt||^2 + (Sg-iQAgAQ - AgAQSg-iQ, AgS/u) 

— ^^qiX^ ‘ ^Il)?Agn) (Ag (daQ'yS^I^Q'yS') 1 Agdl^XLci^ 

— 'y ^ {\Ag\Sg!_^Qa'y\AAQ^fj, AgdpUo) 

k'-9l<5 

y ^ ((Sg' _iQa'-f Sq—lQa'y)AgAgiAQ^fj,AgOi^Ua) 

\q-q\<5 

— y ^ (Ag(Sgi j^ 2 ^Q-yfi^gQa'y)i ^qdfiU(y) (4-4) 

q>q-b 

y ^ {\Ag] S g! A AQ AgdpUa) 

k'-9l<5 

-|- y ^ ((Sg'_^Q'^P Sg—lQ^^^AgA^/AQa-yjAgdj^Ua^ 

\q'-q\<5 

+ E (a.(s,'+ 2A<3 CK7 Ag' Q-yfi), AgdjsUa) + ^(AgdQI AQ); VAgu) 

q'>q-5 

— aX{Ag(\Q\Q),\/Agu) - cX{Ag{\Q\Q tr(Q'^)),VAgu) - ^(AgQ, VAgu) 

12 

“E'^r 

i=i 

Adding up (j4.3p - (l4.4|) and using Lemma lA.ll we have 

Id 10 12 

--(llVAgQIli^ + \\Agu\\l,) + /rllVAgulli^ + rllAAgQIli^ = + E (4-5) 

i=l j=l 

Set 

(fit) := WVQWjjs + ||u|||^s, (flit) := ||S’oVQ||^2 + ||5'on||^2, f2it) := fit) - fiit), 


with fi and f 2 representing the low-frequency part and high-frequency part of f, respec¬ 
tively. Then ()4.5I) leads to the following estimate (see Appendix B for the details): 
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+ r||A,AQ||i,) 

gGN 

<C{l + ||tt||i2||Vu||^2 + ||VQ||^2||A(5||i2 + \\Q\\\2 + ||Q|li4 + \\Q\\\6)‘^{t) 

+ -^W^QWho + ^W'^uWjjs. 

Next, we estimate the low frequencies. Applying Sq to the first equation in (II.2|) . 
multiplying by —SqAQ, taking the trace, and integrating by parts over we have 

— WVSoQWh+rWSoAQWl, 

= {So{u ■ vg), SoAQ) + (SoiQn - QQ), SqAQ) 

- A(5o(|g|hl), SoAQ) + ar(5og, SoAQ) + cr{So{Qtv{Q^)), SoAQ) 

<C\\nh4VQh4SoAQU2 + \\Som-m\\L4So^Q\\L^ 

+ c||5o(|g|^^)||L2||5oAg||i2 + c||5ovg||i2 + c||g||ie||5oAg||i2 
<c||u|||,||Vu|||,||vg|||,||Ag|||,||SoAg||i2 + ||5o(go-og)||ii||5oAg|U2 (4.7) 

+ c||So(|g|f2)||z.i||5oAg||i2 + c||5ovg||i2 + c||g||ie||5oAg||i2 

<c||u||22||Vu|||,||vg|||,||Ag||22||SoAg||^2 + c||g||^2||Vu||^2||5oAg||i2 
+ C\\SoVQ\\12+C\\Q\\14SoAQ\\l2 

<^||5oAg||i2 + c||vg||l,.+c||u||i2||Vn||i2+c||vg||i2||Ag||i2 
+ c||g||i2||vu||i2 + c||g||ie. 

Then we apply So to the second equation in (11.21) . multiply by Sou, and integrate by parts 
over to obtain 

^^ll‘S'o«|li2 +^||5 'oVu|||2 

= -(5o(u • Vu), Sou) + (5o(Vg O VQ), SoVu) - (SoV ■ (QAQ - AQQ),Sou) 

- A(S'oV • (|g|(Ag — aQ - cgtr(g^))). Sou) + ^(SoV • g. Sou) 

< llSo(u-Vn) 11^2 ||5on||i2 + ||5o(vg0vg) 11^2 ||5oVu||i2 
+ |A|||5o(|g|(Ag - ag - cgtr(g2)))||^2||5oVu||i2 

+ ||So(gAg - AQQ)\\l2\\SoVu\\l2 + |n|||SoV • Q\\l2\\Sou\\l2 

< C||So(u • Vu)hi \\Sou\\l2 + c||So(vg O vg)||ii ||5oVu||i2 (4.8) 

+ C||So(|g|(Ag -aQ- cgtr(g2)))||^i ||5oVu||i2 

+ C||So(gAg - AQQ)\44So^uh2 + C\\SoV • Q\\l 2 \\Sou\\l 2 
<C||u||i2||Vu||i2||5ou||i2+C||Vg||i2||5oVu||i2 
+ c(||g||i 2 ||Ag ||^2 + ||g||i 2 + \\Q\\l4\\SoVuh2 + c\\SoVQh 2 \\Souh 2 

< ^||5oVu||i2 + cMls + cwvQWjjs + c\\u\\l, ||Vu||i 2 + c||vg||i 2 
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+ C'(||Q||i2||A(5||^2 + \\Q\\i,2 + ||<5|li4). 

We add ()4.7p - ()4.8p to obtain 

+ ^ll‘S'oVu||^2 + —||S'oA(5||^2 

< Cifit) + C\\u\\l,\\Vu\\l, + C||VQ||i2||AQ||2, + C'||Q||i 2 ||Vn||i 2 (^.9) 

+ C'IIQIIl 2||AQ||^2 + (^(11(511^2 + ||Q||i4 + llQIlie + ||V(5|li2). 

Finally, we derive the estimates of the high norms. Adding (14.6p and (14.9p . we have 

1 ^(() + |||Vu|| 5 ,. + |||AQ|||,. < A{tMt) + B{t), ( 4 . 10 ) 

with 

A{t) =C(l + ||ii|li 2 ||Vu ||^2 + ||V(5||i2||AQ||^2 + ||Q||i 2 + ||Q||i4 + HQUie), 

B{t) =C||u||i 2 ||Vn||i 2 +C||VQ||i 2 ||AQ||i 2 +C||Q||i 2 ||Vu||i 2 

+ C'IIQIlL2||A(5||i2 + C[\\Q\\j^2 + \\Q\\% + \\Q\\% + W'^QWl^)- 

From the a priori estimates in Proposition 12.21 we know that both A{t) and B{t) belong to 
L^(0, T) and increase exponentially in time. Then, by Gronwall’s inequality (Lemma I A. 4p . 
we can conclude that (p{t) increases double exponentially in time as stated in (j4.2p . □ 


5. Weak-Strong Uniqueness in Dimension Two 

In this section, we prove that a global weak solution and a strong one must coincide, 
provided that they have the same initial data {Q,u) G x for s > 0. The 

result is stated in the following theorem. 

Theorem 5.1. Let {Q,u) G x with s > 0, be the initial data. By 

Theorem EU there exists a weak solution {Qi,ui) of system (jl.2p such that 

Qi G L£>,(M+; H^) n lL(M+; H^), ui G LZiR+;L^) n lL(M+; H^). (5.1) 

Theorem KT\ gives the existence of a strong solution {Q 2 ,U 2 ) such that 

Q2 G LS',(M+; 77^+1) nLL(M+; 77^+2), U2 G LS'c(M+; n lL(M+; 77*+'). (5.2) 

Then {Qi,ui) = {Q 2 ,U 2 ). 

Proof. Denote 6Q := Qi — Q 2 and 5u := ui — U 2 . Then {5Q,6u) satisfies the following 
system: 
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'dt 6 Q + 6 u ■ V 6 Q - 6 n 6 Q + 6 Q 6 n + 6 u ■ VQ 2 + U 2 ■ S/ 6 Q + Q 26 n + 5Qn2 - SnQ 2 - n25Q 
= X\Qi\6D + X{\Qi\- \Q2\)D2 
+r(AJ(5 - a5Q - c{ 6 QtT{Ql) + Q 2 tr{Qi 5 Q + SQQ 2 ))), 
dt5u + V{ 6 u ■ V 6 u) 

= fiAdu - V{V • {V 6 Q © V(5Q)) 

< +r{V ■ {dQA 6 Q - AdQ 6 Q)) - V{u 2 ■ V5u + 6 u ■ VU 2 ) 

+P(V • { 5 QAQ 2 + Q 2 A 6 Q - A5QQ2 - AQ 26 Q)) 

-V{V ■ {V 6 Q 0 VQ 2 + VQ 2 0 V^g)) - XV{V ■ (IQiKA^g - a 6 Q))) 

-XV{V • (igii - |g 2 |)(Ag 2 - ag 2 )) + Xcv{v ■ (|gi|<5gtr(g2))) 

+Acp(v • (|gi|g2tr(gi,5g + 5 QQ 2 ))) 

+xcv{v • ((igii - |g 2 |)g 2 tr(gi))) + kv{x 7 • sq). 

(5.3) 

Similarly to the proof of ProDosition l2.ll we multiply the hrst equation in (|5.3p by —A5Q+ 

6 Q, take the trace, and integrate by parts over to obtain 

^^(l|V5g||i2 + ||-5g|li2) - (<5^ • V 6 Q,ASQ) - {5QSn - 6X16Q,A6Q) 

- { 6 u ■ VQ 2 + U 2 ■ XJ5Q + 6 Qn 2 - n 25 Q, A 6 Q) + { 6 u ■ VQ 2 ,5Q) 

- (g2(5fi — 5 QQ 2 , A5Q) + (g2(5n — 5 QQ 2 , sq) 

= -r||A5g||2, - r||V,5g||i. - ar\\V5Q\\l, - ar\\ 6 Q\\l, (5-4) 

+ AdQi |(5iA + (IQi I ~ |g2|).C^2) —A5Q + SQ) + cT[5QtT{Q\), A 6 Q) 

+ cr(g2tr(gi,5g + 5gg2), a5Q) - cr(<5g|gi|2, 6 Q) 

- cT{Q 2 tr{Qi 6 Q + 5QQ2),6Q). 

Multiplying the second equation in (15.3p by 6 u and integrating by parts over yields 

= - (V • {V 6 Q 0 VdQ),du) - {5QA5Q - A 6 QdQ, V 6 u^) 

- {u 2 ■ V 6 u + 6 u ■ Vu 2 , 6 u) — {Q 2 A 6 Q — A 6 QQ 2 ,'S/du^) 

- i 5 QAQ 2 - AQ 26 Q, V5u^) + {V5Q 0 Vg2 + Vg2 0 V 6 Q, V5u^) 

+ X(IQilAdQ,Vdu) - aX(IQil 6 Qi,V 6 u) + X((IQil - lQ 2 l)(AQ 2 - aQ 2 ),V 6 u) 

- cX(lQildQtr(Qf),Vdu) - cA(|gi|g 2 tr(gi(5g + 5 QQ 2 ),X/ 5 u) 

- cA((|gi| - |g2|)g2 - k{6q,v5u). 

Adding (15.41) and (15.51) together, and performing analogous cancellations to those in the 
proof of Proposition 12.11 we have 

^|(liv<5g|ii2 + I|5g|ii2 + \\ 6 u\\i,) + fi\\v 6 u\\i, + rdiA^gni^ + Iiv,5g|ii2) 

= {5u ■ VQ 2 + U 2 ■ V 6 Q, A5Q) + { 5 Qn 2 - ^ 25 Q, A5Q) - {5u ■ VQ 2 , 6 Q) 
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+ X{\Qi\5D, 6 Q) - AdQil - \Q 2 \)D 2 , A 6 Q - 5Q) + cr(5Q tr(Q?), A 6 Q) 

+ cT{Q2tr{Qi6Q + 6QQ2),A6Q) - cT{6Q\Qi\'^,6Q) - cr{Q2ti{Qi6Q + 5QQ2),6Q) 

+ (Su ■ V6u, U2) - {6QAQ2 - AQ26Q, V 5 u^) + (V^g 0 VQ2 + VQ2 O V 5 Q, V6u^) 
-aA(|Qi|5Q,V<5«) + A((|gi| - |g 2 |)(Ag 2 - aga), V^n) - cA(|gi|5gtr(g?), 

- cA(|gi|g2 tr(gi(5g + 5QQ2),v6u) - cA((|gi| - |g2|)g2 tr(g2), vsu) 

- k{5Q,V5u) - ar(||V5g||i. + \\ 6 Q\\l,) 

< \\A 6 Qh. {\\ 6 uh 2 \\VQ 2 \\l^ + \\V 6 Qh 2 \\u 2 \\l^ + 2\\6Q\\^g 
+ PQIIl2 ||vg2||L°° + |A|||(5g||i2 ||gi||Loo||v5tt||2,2 

+ |A|(||A<5g||^2 + ||5g||^2)pg||^ll|i)2||^^+cr||A5g||i2||5g||i4||gi||23 
+ cr||A,5g||i2||g2||L-(||gi||L4 + ||g2||L4)||<5g|k4 + cr||<5g||i4||gi||L2||gi|Ucx> 

+ cr||(5g|||4||g2||L2(||gi||Loo + ||g2||L°°) + II^^2 ||l°°||(^w||j;^2||V(5m||£2 

+ 2||<5g|| 2 ||Ag 2 || 2 i|v5u||i2+2||v5g||i2||v<5n||i2||vg2||L- 

L/ s l — s 

+ |aA|||V5u||^2pg||^2||gi||i^ + ci|V5«||i2(||5g||^ll|Ag2||^^ + ||5g||i2||g2||L^) 

+ c|A|||V5tti|j;^2 ||(5g||j;^4 (||gi||£,oo(||gi|||8 + ||g2|||8) + ||g2|li8||g2||L°°) 

+ \k\\\V 6 u\\l 2 \\ 6 Q\\l 2 - ar(||v<5g||i. + ||<5g||i.) 

<^\\A 6 Q\\l, + i^\\V 6 u\\l, + C(||Vg2||ioo + \\u 2 \\l^)\\Su\\l, 

+ ^(1 + ||g2||ioo + iigiiiioo) \\ 6 q\\ 1 , + c(i + ||u2||ioo + ||vg2||ioo)||v5g|ii. 

+ c'r(iigiii^8 + iig2iii°o(iigiii|4 + Iig2iii4) + iigiiiisiigiiiioo + iig2iii8iig2iiioo 

+ iiQiiii°°(iiQiiiL8 + iiQ2iii8) + (iigiiiL°° + iig2iiL°°)(iigiiiL2 + iig2iiL2)) ii<^giii4 

+ C{\\Vu 2 f 2 +\\AQ 2 f 2 )\\SQ\\\_. 

Since we restrict ourselves to the two-dimensional case, ||5g|||4 and ||5g|p2 can be con- 

L s 

trolled by \\5Q\\l, + \\V 6 Q\\l,. Moreover, we know from the imbedding theorem that 

1|V'U2|| 2 < C'||u2||h'>+i, ||Ag 2 l| 2 < C\\Q 2 \\h‘+ 2 - 

Li-s L i-s 

In addition, from the conditions of gi,tti,g 2 , and U 2 , he., (I5.ip - (l5.2p . we know that all 

the coefficients of ||(5 m||? 2, ||<5g||?2) ||V(ig||? 2 j ||^g||? 4 ) and ||(ig|p 2 are integrable with 

L s 

respect to time. Therefore, we can use Gronwall’s inequality (Lemma IA.4I) to conclude 
the uniqueness of the solution. □ 


Appendix A. Some Basic Theories and Lemmas 

In this section, we review some important theories and lemmas that are used extensively 
in this paper. 

First, let us introduce the Littlewood-Paley theory; see [2] for the details. 
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Proposition A.l (Dyadic Partition of Unity). Let C be annulus G : | < |^| < |}. 
There exist radial functions x o.n-d valued in interval [0,1], belonging to |)) and 

T){C), respectively, such that 

Xii) + X] = 1 for any i G M"*, 

i>i 

= 1 for any ^ €W^\{0}, 

jez 

\j - j'l > 2 ^ supp{(f{2~^■)) n supp(<y9(2"^ •)) = 0, 
j >l ^ supp(y) n supp((/?(2"^-)) = 0. 

Definition A.l. The homogeneous dyadic blocks Aj and the homogeneous low-frequency 
cut-off operators Sj are dehned for all j G Z by 


AjU = jr ^(vj (2 ^f)Fu) = 2 -^'^ f h{ 2 ^y)u{x — y)dy, 

SjU = J='-^{x{ 2 ~^f)J^u) = 2 ^"* [ h{ 2 ^y)u{x - y)dy, 

where J- denotes the Fourier transform on and 

h := h := T~^X- 

We define the Sobolev norm of space as 

:= {\\Sou\\l2 + '^2^‘^^\\Aju\\l2)^. 

jeN 

Next, we recall the Bony’s paraproduct decomposition for two appropriately smooth 
functions u and v: 

uv = TuV -|- TyU + R{u, v) = TyV + T^u, 


TuV Sj-iuAjV, R{u,v) := AjuA-iv, 

J b-/l<i 

TyU = TyU + R{u, v) = Sjj^2vAjU. 

j 


where 
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Then we have 

Aj{uv) = AjTuV + AjT^u 

= Sy_^uA.>v + Aj S.^^vA.u 

3 3 

= Y ^ 3 (Sj'_iUA.,v)+ Y 

b-i'l<5 3 > 3 -^ 

= Sj-iuAjV+ Y^ _^u]A.iv 

Ii-i'l<5 i'>i-5 


In order to obtain the uniqueness and higher regularity of the weak-strong solutions in 
§4-§5, we now recall the following useful inequalities that we use extensively. These in¬ 
equalities follows from the Bernstein-type lemma and the commutator estimates in Chapter 
2 of [2] by the construction of Aj. 


Bernstein-type inequalities: 

IIVS’juIIlp < C'2'^ IIuIIlp for any 1 < p < oo, 

C2^ WAjuWiP < ||VAju||LP < (72-^ II AjulliP for any 1 < p < oo, 

IIAjuIIl? < \\Aju\\Lv with 1 < p < g, 

Il'SjUtlli? < \\Sju\\LP with 1 < p < g. 


Commutator estimates: 

||[Aj,u]u||i;,r < C'2 "'^'||Vu||lp||u||l9 


■ . 1 1 

with —I— = 
p q 


1 

r ’ 


where C is independent of p, q, and r. 

In order to deal with the highest derivatives of u in the hrst equation of system (11.21) 
(and the corresponding one in the approximation systems) and the highest derivatives of 
Q in the second equation of the system (and the corresponding one in the approximation 
systems), we need to introduce the following lemma: 


Lemma A.l. Let Q and Q' be two d x d symmetric matrices, and let Ll = ^(Vu — Vtt^) 
be the vorticity with {S/u)ai 3 = dpUa- Then 

[LlQ' - Q'n,AQ) - (V • {Q'AQ - AQQ'),u) = 0. 


Proof. For any two d x d matrices A and B, we know that tr{AB) = tT{BA). Then we 
have 

{LIQ' - Q'n, AQ) - (V • {Q'AQ - AQQ'),u) 

= (LIQ' - Q'n, AQ) + {Q'AQ - AQQ',Vu^) 

= inQ',AQ) - {Q'n, AQ) + {Q'AQ - AQQ',Vu^) 

= {Q'AQ, n) - {AQQ',n) + {Q'AQ - AQQ',Vu^) 
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= {Q'AQ-AQQ',n + Vu^) 
= iQ'AQ-AQQ',D) = 0, 


where, in the last equality, we use the fact that Q, Q' , and D are symmetric. □ 

Remark A.l. For any d x d matrix Q, we have 

tr(Q^) < ^|tr((5^)P + ^tr((5^) for any e > 0. (A.l) 

Proof. Let x, y, and 2 ; be the eigenvalues of Q. Then we have 

tr{Q)=x + y + z, tr((5^) = tT{Q^) = x^ + y^ + . 

Since it is obvious that tr(Q'^) < |tr((5^)p, we only need to show that, for any e > 0, the 
following inequality is true: 

tr(Q^) < ^tr(Q^) + -tr(Q2). 

o £ 


By Young’s inequality with e, we have 

„2 


q X £ A 

x^ < -h -X^ 

~ £ 4 




Then we obtain the desired inequality. 


□ 


Lemma A.2 (Aubin-Lions lemma). Let Xo,X, and Xi be three Banach spaces with Xq C 
X C Xi, let Xq be compactly embedded in X, and let X be continuously embedded in Xi. 
For 1 < p,q < 00 , let 


W= {u€ LP(0,T; Ao) : u G L''(0,T; Ai)}. 


Then 

(i) If p < 00 , then the embedding of W into L^(0,T; Ai) is compact. 


(ii) If p = 00 and q > 1, then the embedding of W into (7(0, T; A) is compact. 

Lemma A.3 (Gagliardo-Nirenberg interpolation inequality [l3]). Let 1 < q,r < 00 . For 
0 < j < m, the following inequalities hold: 


where 


1 j /I n1 

- — —h a(-—) + (1 — a)-, 

p n r d q 


for all a in the interval 

— < a < 1, 
m 

and C = {d,m,j,q,r,a), with the following exceptional cases: 

(i) If j = 0, rm < d, and q = 00 , then we make the additional assumption that either 
u tends to zero at infinity or u a L^ for some finite s > 0. 

(ii) If 1 < r < 00 , and m — j — I is a nonnegative integer, then it is necessary to 
assume additionally that a ^ 1. 
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Lemma A.4 (Gronwall’s inequality). Let a > 0 and fd > 0 be integrable functions on 
[0,T]. If a differentiable function Y satisfies the differential inequality: 

Y'{f) < a{t)Y{f) + fl{t) for t € [0,T], 

then 

y{t) < 1^(0) exp ^ J a{s)ds'^ + j /3(s)exp ^ J a{T)dT^ ds for any t € [0,T]. 


Appendix B. The Estimates of Inequality (14.61) 

Before proving inequality (|4.6p . let us first introduce the following Lemma. 

Lemma B.l. Let u G D L^ with p > 1 and s > 0. Then, for any k >2 and g G N, 






(B.l) 


where {aq^k{t)}q&fi is a sequence in P. 

Proof. We prove this lemma by induction. 

Firstly, for A: = 2, by using the Bony’s paraproduct decomposition, we have 

Aq { u ^) = Sq-luAqU + ^ [Ag, 5^/U 

l<?-g'l<5 

+ Y1 5’g_iu)AqA^/u+ Aq{Sgr^.^uA 

\q-q\<5 q'>q-5 

= E 


.'U 


l<i<4 


Let us calculate the right side term by term as follows, 

W^iWlp = \\Sq-luAqU\\Lv < C\\u\\Lv\\AqU\\L^ < C\\u\\Lp\\AqY u\\l2 
<C2-^^4\t)\\u\\Lp\\Vu\\Hs, 
where {aq^^(t)}qgp} is a sequence in P, 

||2i2||LP = II [Aq,5g/_^tt]Ag/tt||i:,p < C Y^ 2“^||V5^/_j^m||lp||A^/-u||lo 

\ q -< i '\<5 k-g'l<5 

<C Y 2'''“''lkl|Lp||A^,V«||i2 

l<?-g'l<5 

<C2-i^ Y 2(^'-‘')(^-")a(!^(t)||u||Lp||Vn||i^. 
k-<?'l<5 

where {df\t)}q(z^ = {^q-q (^)}gGN is a sequence in P, 

II^sIIlp = II Y - ‘S'q-lu)AqA^/u| 

I<?-9'|<5 


\LP 
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<C ^ \\S^>_^U- Sq-iu\\Lp\WA^fu\\L^ 

\q-q'\<5 

k-<?'l<5 

\q-q\<5 

<C2-'J^aP{t)\\u\\Lp\\Vu\\Hs, 

with {4^^(t)}g6H = ^''a^^\t)}qen e 

II^4||l. = || E Ag('S'q'+2^A^'t()||LP < C E l|5,'«'‘ll«IIE«lk» 

q'>q—5 q'>q—5 

<C Yl \\u\\lp\\A^^Vu\\l 2<C2-‘^^ Y ^''‘^-^'^^\\u\\Lpa^]\t)\\Vu\\Hs 

q'>q—5 q'>q—5 

<C'2-^^af (t)||n|Up||Vn||H=, 

with {a‘^^\t)}q^^ = {Eg'>g -5 € P. 

Then taking aq^ 2 {t) = maxi<fc< 4 {aq^^(t)}, we know that the result holds for k = 2. 
Next, we first assume the statement is valid for k and then check the validity for the 
case k + 1. Similarly to the previous steps, we have 

Aq{u^+^) = Aqiu^u) = Sq.lU^AqU + ^ [A„ Sg^_Y]Aq'U 

\q-q'\<5 

+ Y iSq'_^u’^-Sq.^U^)AqA^,U+ Y 

\q-q'\<5 q'>q-5 

= E 

l<i<4 

and 

WMlp = WSq-lu’^AqUUp < C\\u’^\\Lp\\AqU\\Loo < C11U || ^,, || A, Vu 11 ^2 

\\J2\\lp = \\ Y [Aq,S^,_Y]Aq'U\\LV<C Y ‘^-n^S^,_^v!^\\Lv\\A^'u\\L^ 

\q-q\<^ k-<?'l<5 

<C 2^-^u^\\lp\\A^,Vu\\l2 

k-<?'l<5 

<C2-’* E 2'’'-’>'‘-'>aW(()ll«lltJ|Vu||f,. 

k-g'l<5 

<C2-^^af{t)\\u\\l,,\\Vu\\Hs, 
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l<?-g'l<5 

<C Y. WSg'.^v!^ - s^.^u^WlAWAl^ 

k-<?'l<5 

<C Y \\A\LA\^q^u\A 

k-9'l<5 

<C2-<i^-A^\t)\\u\\Y\\^u\\Hs, 

and from the induction assumption, one has 

||J4||Lr = || \{S ,^^uI^ ,uA\\l.<C Y hhAAg'uA ^ 

^ .. . L^-1 



q'>q—5 


q'>q—5 

< 

c Y 

\ikp2 ^ ^0,^'^^{t)\\u\ 

\Tkp\MH^ 


q'>q—5 



< 

C2-'?* Y 


iUw^Ahs 


q'>q- 

5 


< 

C2-'?*ag,fc+] 

Xt)\\u\\lA\^u\\Hs, 



with {ag^k+l{t)}qm = {Eg'>g -5 2^'^ ^ ^ 

Then taking aq^k+i{t) = max{aq^^(t), (t), Oq^fe+i(t)}, we know that the statement is 
true for the case k + 1. By induction, we complete the proof. □ 

Then, in order to derive (14.6h . we estimate the terms on the right-hand side of (j4.5|) 
one by one. For the detailed estimates of the following terms, we refer the readers to the 
appendix in [47] : 

\IA < C2-25^6g(t)(||n|||,||Vu|||,||VQ||i||AQ|||. 

+ ||VQ|||,||AQ|||,||u||i||Vu||i||Ag||H^), 

IX2I + IX3I + IX5I + iXel < C' 2 - 2 ''* 6 g(t)||VQ||| 2 ||AQ||| 2 ||n|||.||Vn|||.||AQ||^^., 

IX 4 I + \IA < C2-2'?^6g(t)||u|||,||Vn|||,||VQ|||.||AQ|||., 

|:7i| < C2-29^6g(t)||u|||,||Vu|||,||u||i||Vu||i, 

1^21 < C2-2'?^6g(t)||VQ|||,||Ag|||,||Q||i||Vg||i||Vu||H=, 

IJ 3 I + \JA + \JA + m < C2-29^6g(t)||Vg|||,||AQ|||,||u||i||Vn||i||AQ||^^., 

\JA + \Js\ < C2-2''*6g(t)||VQ|||,||AQ|||,||VQ||i||Ag|||.||Vu||H», 

where {bg{t)}g^n is a sequence in A. 

Next we will give the details of the estimates for the remaining terms: 

IX 9 I = |ra(AgQ,AAgQ)| =ra||AgVQ||2, <C2-2^^6W(t)||VQ|||,., 

IJ 12 I = |K(AgQ, AgVn)| = |K(AgVg, Agu)! 
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< C\\AgVQ\\L2\\A,u\\L2 < C2-^<^%f\t)\\VQ\\Hs\\u\\Hs, 

for some {bq\t)}q^^ S with i = 1,2. Applying Lemma IB. 11 we have the following 
estimates: 

iXiol = |rc(A,(Qtr(g2)), AA,g)| < C\\Aq{QtT{Q^))\\L2\\AAqQ\\L2 
<C2-'?X3WIIQIIi4||vg||j^.2-'^^aW(t)||Ag||H^ 
<C2-^‘^^bf{t)\\Q\\l4VQ\\H4m\Hs, 

\Jw\ = \ - aX{Aq{\Q\Q),VAqu)\ 

<C\\Aq{\Q\Q)\\L2\\\/AqU\\L2 

<C2-'?X2WIIQIlL2||Vg||H»2-'?^aW(t)||Vu||H» 

<Cis)2-^^%^\t)\\Qh2\\VQ\\H^\\\/u\\H^, 

\Jn\ = I - cA(A,(|g|gtr(g)2), VA,n)| 

<c||A,(|g|gtr(g2))||^2l|VA,u||i2 

<C'2-^X4(i)||g|li6||Vg||H»2-''^a«(t)||V«||H» 
<C2-^<i^bf{t)\\Q\\l4VQ\\H4^u\\Hs, 
for some {6q*^(t)}ggp} G with i = 3,4,5. 

Finally, for the remaining two terms Xg and Xg it is not easy to obtain the bounds 
directly. The key idea is to combine them together and identify the cancellation factors 
between them for the estimates. Applying the Bony’s decomposition, we have 

X8 = -A(A,(|g|ii),AA,g) 

= -A(5,_i(|g|)A,AAA,g)-A ([A,;5^,(|g|)]A^,AAA,g) 

k'-g|<5 

-A ((5^._jg|-5,_i|g|)A,A^,AAA,g) 

l<?'-g|<5 

-A ^ (A,(5^,^2^A^'IQI),AA,g) 

q'>g-5 

4 

= A^Xg^j, 
i=l 

Jg = X{Aq{\Q\AQ),VAqU) = X{Sq-l\Q\AqAQ,VAqU) + X{[Aq-,S^>\Q\]A^>AQ,VAqU) 

+ A ^ {{S^^_,\Q\ - Sq.i\Q\)AqA^,AQ,VAqu) 
k'-g|<5 

+ A iAqiS^,^^AQA^.\Q\),VAqu) 

q'>q-5 

4 

= A^Xg^i. 
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Since Q is symmetric and = Vir — D is skew-symmetric, we find that, for 1 < f < 3, 

^8,i + <J9,i = 0- 

Then 

Tg -I- Jg = Is,A + J9,4- 

In addition, we have 

|X8,4| = | ^ {A,{S^,^,DA^,\Q\),AAgQ)\ 

q'>q-5 

E \\Sg'+2^u\\L4\'\Q\\\L4^^qQ\\L^ 

q'>q-5 

<c 2'^'+2||^||^,2-'''||A^,VQ||i4||AA,Q||^2 

9 >9-5 

E lkllil|Vrx|||,||A^,VQ|||,||A^,AQ|| 2 ,||A,AQ|| 

9 >9-5 

< _ 

9 '>9-5 


HO 


(, 2 - 2 .. ^ 2<»-'>')*s<;)(t)4i){()||„||j,||v«|||.||VQ||j.||AQ||i. 


< C2-"’*6l«l (() Hull LII Vu|| I IIVQII ’ IIAQII ’ , 


and 


1X9,41 = 1 E i^q{S^'+ 2 ^Q\\Q\),'^^q 4 \ 

9'>9-5 

<c E i>, I _^2^Q\\l'>‘ 11‘^q' QIIl4 

9 '>9-5 

<C 2'''||VQ||i42-^'||A^2VQ||i4||A,Vu||,.2 

9 '> 9-5 

<C Y l|VQ||il|AQ|||,||A^2VQ||E|A^2AQ|||,||A,Vn||^2 

9 '> 9-5 

< C ^ \\VQ\\l4AQ\\l2-'i'^a^^]\t)\\VQ\\^^^^^ 

q'>q-5 

<C2-^qsbf\t)\\VQ\\l4AQ\\l\\VQ\^^^^^^ 

with {bf\t)}q(.n = {Eg'> 9 - 52 ^''"'' ^''a^^\t)aq\t)}q(.n G l^. 

Multiplying all the above estimates by adding them together, taking the sum in q 
for g G N, noticing that {6g(t)}ggN, {bq\t)}q^f^ G with 1 < f < 6 , and using the Cauchy 
inequality with suitable e, i.e., ab < ea^ -|- ^b"^, we obtain the desired estimate (j4.6l) . 
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